#ifndef _UTILS_
#define _UTILS_

#include <myTypes.h>
#include <myMacros.h>
#include <assert.h>

double myMse (dVec_t a, dVec_t b) {
  double res = 0;
  assert (a.size() == b.size());
  forEach (i, 0, a.size())
    res += (a[i] - b[i]) * (a[i] - b[i]);
  return res / a.size();
}

#endif
